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account  for  projectile  mass  loss  effects.  The  current  model  does  demonstrate 
good  correlation  with  ballistic  test  results  in  the  standard  velocity  range. 
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Abstract 


The  research  program  consisted  of  two  main  parts:  (1)  a  series  of 
ballistic  tests  with  long  rod  penetrators  of  SAE  4140  steel  impacting 
target  plates  of  RHA  steel  and  5083  aluminum  alloy,  and  (2)  the 
development  of  an  improved  five  stage  engineering  model  of  the  ballistic 
perforation  process.  In  addition,  a  spl it-Hopkinson  (Kolsky  bar) 
arrangement  was  set  up  to  measure  material  strength  at  high  rates  of 
loading.  Results  of  the  present  ballistic  tests  will  be  used  in 
conjunction  with  subsequent  improvement  of  the  engineering  model  to 
account  for  projectile  mass  loss  effects.  The  current  model  does 
demonstrate  good  correlation  with  ballistic  test  results  in  the 
standard  velocity  range. 
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Description  of  Research  Program 

The  overall  research  program  consisted  of  two  major  parts  and  a 
secondary  supporting  project.  The  first  major  task  was  the  design 
and  performance  of  a  series  of  ballistic  tests  with  long  rod 
penetrators  made  of  SAE  4140  steel.  These  were  used  to  impact  target 
plates  of  rolled  homogeneous  armor  steel  (RHA)  and  5083  H- 115  aluminum 
alloy.  Thicknesses  of  the  RHA  plates  were  75.0  mm  and  50.8  mm  and  the 
aluminum  target  consisted  of  3  plates  in  contact,  each  of  44.5  mm,  for 
a  total  thickness  of  133.5  mm.  Impact  was  normal  to  the  plate  with 
impact  velocities  ranging  from  1356  to  1560  m/s.  Eight  of  the  shots 
were  successful  in  supplying  useful  information. 

The  tests  were  conducted  in  a  well  instrumented  ballistic  range 
which  provided  pre-impact  velocity  measurements  by  photo-cells  and 
X-ray  photographs  of  the  projectile  immediately  prior  and  after 
perforation  by  flash  radiographic  equipment.  After  completion  of  the 
ballistic  tests,  the  target  plates  were  sectioned  and  etched  around 
the  perforation  cavities  which  were  then  subjected  to  close  inspection. 
Mechanical  property  tests  were  also  made  on  the  projectile  and  target 
materials.  Since  the  perforation  process  occurs  at  very  high  rates  of 
loading,  mechanical  properties  at  those  rates  would  be  useful  in  the 
analytical  modelling.  This  data  can  be  obtained  by  means  of  a  split 
Hopkinson  pressure  bar  (Kolsky  bar). 

When  the  measured  residual  velocities  following  perforation  are 
extrapolated  to  zero  by  an  analytical  procedure,  an  estimate  could  be 
obtained  for  the  ballistic  limit  velocity,  VL>  for  the  projectile  and 
the  various  target  plates.  For  the  steel  plates,  for  the  75.0  mm 
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thickness  was  obtained  to  be  1325  m/s,  while  that  of  the  50.8  ran  steel 
target  plate  was  V|_  =  1150  m/s.  In  the  case  of  the  aluminum  targets, 
the  test  data  was  not  sufficient  for  a  reasonable  estimate  of  V^. 

Some  of  the  more  important  results  of  the  various  tests  and 
examinations  are  attached  to  this  report  as  Appendix  A.  The  attached 
figures  and  tables  consist  of  the  following: 

Fig.  1  -  Geometry  of  projectile  of  SAE  4140  steel. 

Table  I  -  Penetrator:  Material  composition. 

Table  II  -  Penetrator:  Mechanical  properties. 

Table  III  -  Target  plates:  Mechanical  properties. 

Table  IV  -  Projectile  velocity  characteristics  prior  to  impact. 

Table  V  -  Contact  time,  residual  velocity,  and  physical  and 
geometrical  characteristics  after  perforation. 

Table  VI  -  Dimensions  of  projectile  entrance  and  exit  holes  in 
the  target  plates. 

Figs.  2,3  -  Cavity  profiles  in  steel  target  plates  showing 
hardness  (Rc)  measurements. 

Fig.  4  -  Cavity  profile  in  an  aluminum  target. 

Fig.  5  -  Hole  profile  after  perforation  of  50.8  mm  RHA 

target. 

Fig.  6  -  Hole  profile  after  perforation  of  75  mm  RHA  target. 

Fig.  7  -  Hole  profile  after  perforation  of  133.5  mm  AL  5083 

H- 115  target. 

Fig.  8  -  Hole  profile  after  perforation  of  133.5  mm  AL  5083 

H-115  target. 

Fig.  9  -  Characterization  of  the  final  deformation  and  hole 

dimension  for  Table  VII. 

Table  VII  -  Characteristic  parameters  measured  for  defining 
the  final  target  deformation. 

This  examination  of  the  results  of  the  ballistic  tests  is 
intended  to  serve  as  a  physical  basis  and  motivation  for  the  development 
of  an  engineering  model  of  tlfe  perforation  process  at  high  impact 
velocities.  The  test  results,  especially  the  residual  velocities  and 
residual  projectile  lengths,  will  also  serve  as  reference  data  for 
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subsequent  comparison  with  analytical  predictions.  The  model  would 
include  the  effects  of  erosion  of  the  projectile  and  of  plastic 
deformation  and  inertia  of  the  projectile  and  target. 

The  second  major  task  was  the  development  of  a  new  five  stage 
engineering  model  of  perforation  of  a  metal  target  by  a  rigid 
projectile  of  standard  ordnance  velocities.  The  initial  stage  relies 
on  an  assumed  velocity  field  for  a  plastically  deforming  medium  with 
a  moving  boundary  (the  indenter),  where  the  unknown  parameters  are 
determined  from  the  upper  bound  theorem  of  plasticity  modified  to 
include  dynamic  effects.  Subsequent  deformation  stages  include  the 
development  of  a  bulge  at  the  rear  surface  of  the  target,  the  extension 
of  the  bulge,  the  development  of  a  plug  and  a  shear  zone,  and  the 
ejection  of  the  projectile  and  plug  following  material  failure.  This 
model  is  an  improvement  over  that  of  Awerbuch  and  Bodner  as  it  does 
not  rely  on  any  empirical  factors  and  is  more  suitable  for  extension 
to  the  case  of  projectile  erosion.  A  paper  describing  this  model  has 
been  prepared,  "Dynamic  Perforation  of  Viscoplastic  Plates  by  Rigid 
Projectiles",  by  Ravid  and  Bodner,  and  appears  as  Appendix  B  of  this 
report.  It  has  been  accepted  for  publication  in  the  International 
Journal  of  Engineering  Science  (for  publication  in  1983).  Further 
improvement  of  the  model  to  include  the  effects  of  projectile  mass 
loss  will  lead  to  a  theory  that  could  be  correlated  with  the  results 
of  the  test  program. 

The  secondary  supporting  project  was  the  setting  up  of  a  split- 
Hopkinson  bar  arrangement  (Kolsky  bar)  in  order  to  obtain  stress- 
strain  relations  at  high  rates  of  loading.  The  experimental  facility 


is  shown  schematically  in  Fig.  1.  A  12.7  mm  diameter  steel  projectile 
embedded  in  a  teflon  collar  is  propelled  in  a  smooth  bore  (18  mm) 
powder  gun  and  strikes  the  input  bar  which  is  made  of  the  same  material 

as  the  projectile  (PH  -  15-5).  Both  input  and  output  bars  are  identical, 

each  12.7  mm  in  diameter  and  876  mm  long.  The  length  of  the  projectile 
in  the  present  tests  was  100  mm  and  its  impact  velocity  was  in  the  range 
of  15  to  35  meters  per  second.  The  maximum  stress  produced  by  the  higher 
velocity  was  about  720  MPa,  which  was  a  little  over  half  the  yield  stress 
of  the  Hopkinson  bars.  Such  low  velocities  were  obtained  by  shortening 
the  cartridges  and  using  between  0.5  and  1.0  grams  of  gun  powder.  The 

specimens  were  short  disks,  12  mm  in  diameter  with  aspect  ratios  U/d)  of 

0.3  to  0.5.  The  strain  pulses  were  detected  by  two  pairs  of  strain  gauges 
(WK-06-125AD-350  by  Micro  Measurement)  and  recorded  by  a  two  channel 
digital  oscilloscope  (Explorer  204  by  Nicolet). 

Four  metals  and  four  polymers  were  each  tested  at  two  or  three  impact 
velocities  and  the  tests  were  repeated  at  least  three  times.  The  materials 
tested  were: 

(1)  SAE  1070  Steel 

(2)  2011-T3  Aluminum 

(3)  99.93 %  Purity  Electrolytic  Copper 

(4)  Commercial  Brass 

(5)  Polyamid  (Akulon)  6 

(6)  Polytetraf luoroethylene  (Teflon) 

(7)  Structural  adhesive  FM-300K 

(8)  Structural  adhesive  FM-73 

The  oscilloscope  records  of  the  tests  conducted  with  the  non  metallic 
specimens  showed  that  the  reflected  pulses  were  large  compared  to  the 
transmitted  ones,  and  since  the  strain  rate  is  approximately  proportional 
to  the  magnitude  of  the  reflected  pulse,  this  means  that  the  strain  rates 
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in  the  tests  of  the  non  metallic  specimens  were  higher  than  in  the 
metallic  specimens  (where  the  reflected  pulses  were  relatively 
smaller).  It  was  also  clear  from  the  oscilloscope  records  that  the 
strain  rates  were  not  constant  throughout  the  tests  but  varied  in 
time.  In  the  tests  of  the  non  metallic  materials  one  could  relate 
to  average  values  of  strain  rates,  whereas  in  the  tests  of  the 
metallic  specimens  average  values  were  meaningful  only  over  small 
portions  of  the  loading  pulses. 

Some  of  the  stress-strain  curves  are  shown  in  Figs. 2-5.  The 
curves  of  the  non  metallic  specimens  are  given  with  indications  of 
the  average  strain  rates  for  each  curve,  whereas  the  results  of  the 
metal! ic  specimens  are  presented  with  indications  of  the  maximum 
values  of  the  strain  rates. 

The  personnel  who  worked  on  the  research  program  during  this 
period  are:  Professor  S.  R.  Bodner,  Associate  Professor  J.  M.  Lifshitz, 
and  Mr.  Moshe  Ravid  and  Mr.  Ezra  Scher,  both  graduate  students. 


CH.1 


Schematic  diagram  of  a  split  Hopkinson  bar  arrangement  and  instrumentations. 
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Typical  results  of  stress-strain  curves  from  "static"  tests  and  under  high  rates  of 
strain  -  FM-73. 
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•'  Typical  results  of  stress-strain  curves  from  "static"  tests  and  under  high  rates  of 
strain  -  Aluminum  2011-T3. 


Table  I  :  Penetrator  Material  Description.  SAE  4140  [EN  19A] 


Symbol 

C 

Si 

Mn 

D 

S 

D 

Mo 

[%1 

BS 

n.2l 

n.72 

0.015 

CNJ 

CN 

C 

c 

1 .06 

n.22 

Table  II  :  Mechanical  Properties  of  Penetrator  Material. 


Ultimate  Stress 

E longat i on 

Impact  Test 

U.T.S. 

e 

•  ZOD 

SW99R 

max 

■M 

■DM 

% 

„  f 

Kg  -m 

78  -  8o 

90  -  91 

22  -  24 

11  -  13 

Table  III:  Target  Plates  and  Mechanical  Properties  of  Target  Material. 


Plate 

Material 

R  .  H  A  .  ^  ' 

Plate 

Th i cknes  s 

(*) 

Mechanical  Properties  - 

Yield  Stress 

Ultimate 

Stress 

Max  . 

Strain 

T 

VV2^ 

U.T.S. 

! 

'max 

mm 

mssmm 

Kg^/mm^ 

50.8 

86 

97-98.7 

15 . 1 -15.3 

(*■'() 

R.H.A. 

75 

80-85 

91 .3-98.6 

15.3 

Ai-5083  H-115 

At,. 5 

1  t-12 

26.2-28.2 

15.8-18.3 

(*)  At  Strain-Rate  of  n-O.I  Sec  '  for  Steel  and  "v 0.6  Sec  '  for  Aluminum 


Plates(obtained  on  M.T.S.  machine,  25  ton  capacity). 

(**)  According  to  Mil-S-12560B  (ORD) 


Table  !V  :  Summary  o'  Before  Tarle;  Penetration  Results  and  Analysis 
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A15083  1 1391  1387  1388  I  -0.61*  0.81*  .1  0.12 
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Maximum  diameter  of  plastically  deformed  forward  end  of  projectile. 


Table  VI  .  Summaiy  of  target  Hole  Oirienb>cn  Heasurement j 


Maximum  diameter  of  circumferential  fracture  region  surrounding  exit  hole  (circumferential  avg.). 
Depth  from  rear  surface  of  circumferential  fracture  region  surrounding  exit  hole  (circum.  avg.)- 
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b.  Enlargement  of  exit  region  (x2.5)  c.  Enlargement  of  entrance  region  (x2.5) 


Etching:  Boiling  50%  HC1 

Figure  5:  Hole  profile  after  perforation  of  50.8  mm  R.H.A.  target. 

(Round  No.  15;  V  =  1461  ro/sec;  V„  =  1310  m/sec;  LD  =  129  mm) 

S  K  K 
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c.  Enlargement  of  entrance  region  (x2.5)  b.  Enlargement  of  exit  region  (x2,5) 


Etching:  Feric  Chloride  +  HC1  +  Alcohol 


Figure  8:  Hole  profile  after  perforation  of  133.5  mm  At,  5083  H-115  target. 

(Round  No.  19;  Vs  =  1560  m/sec;  VR  =  1390  m/sec;  Lr  =  138  mm). 


DYNAMIC  PERFORATION  OF  VISCOPLASTIC  PLATES 
BY  RIGID  PROJECTILES* 

by 

M.  Ravid  and  S.  R.  Bodner 
Faculty  of  Mechanical  Engineering 
Technion  -  Israel  Institute  of  Technology 
Haifa  32000  Israel 


Abstract 


A  model  is  formulated  for  the  dynamic  perforation  of 
viscoplastic  plates  by  rigid  projectiles.  The  process  is 
considered  to  occur  in  five  continuously  coupled  but  distinct 
stages  which  are  amenable  to  analytical  treatment.  An  essential 
feature  of  the  analysis  is  the  use  of  postulated,  physically 
motivated,  deformation  mechanisms  in  conjunction  with  the  upper 
bound  theorem  of  plasticity  theory  which  is  modified  to  include 
dynamic  effects.  Special  attention  is  given  to  the  bulging 
process  and  effects  associated  with  the  later  stages.  This 
model  is  self-contained  and  capable  of  predicting  the  exit 
velocities  of  the  projectile  and  the  plug.  It  also  determines  the 
bulge  and  plug  shape,  provides  the  force-time  history  of 
the  process,  and  describes  a  number  of  geometrical  features 
of  the  transient  and  final  deformation  state  of  the  target 
plate . 


*The  research  work  reported  in  this  paper  has  been  sponsored 
in  part  by  the  U.S.  Army  (European  Research  Office,  London) 
under  Contract  DAJA37-81-C-0047. 
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Introduction 


During  the  penetration  and  perforation  process  of  a 
rigid  cylindrical  rod  in  a  metallic  target  plate,  the  target 
material  undergoes  extensive  plastic  deformation  due  to  the 
high  pressure  developed  by  the  penetrator/ target  interaction. 
After  a  short,  high  pressure  impact  phase  which  could  exist  for 
about  a  microsecond,  a  stable,  steadily  growing  plastic  flow 
field  develops  in  the  target.  This  flow  field  has  a  major 
influence  on  the  penetration  process  since  it  governs  the 
resistance  to  penetration  in  the  initial  stage  prior  to  the 
onset  of  bulging  and  is  also  an  important  factor  on  the  exit 
of  the  projectile.  In  previous  engineering  models  of  pen¬ 
etration,  e.g.  [1],  the  static  and  dynamic  effects  of  the 
flow  field  have  been  considered  in  a  very  approximate  and 
overall  manner  which  required  the  use  of  some  empirical 
factors.  In  addition,  the  dynamic  effects  that  have  to  be 
taken  into  account  in  an  overall  approach  are  not  entirely 
obvious  as  has  been  discussed  by  Tate  [2] .  Another  limitation 
of  the  simplified  model  is  the  difficulty  of  properly  incor¬ 
porating  the  influence  of  material  strain-rate  dependence. 
Despite  these  shortcomings,  the  simplified  three  stage  model 
of  [1]  has  been  reasonably  successful  in  predicting  the  most 
important  features  of  the  perforation  process  [7]  . 
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In  this  paper,  the  process  is  analyzed  by  means  of  a 
five  stage  sequence  of  deformation  mechanisms  shown  in  Fig.l. 
These  are  1.  dynamic  plastic  penetration,  2.  bulge  formation, 
3.  bulge  advancement,  4.  plug  formation  and  exit,  and 
5.  projectile  exit.  The  first,  principal,  stage  involves 
the  determination  of  an  approximate  plastic  flow  (velocity) 
field  surrounding  the  projectile  at  each  increment  of 
penetration  and  includes  the  development  of  a  "lip"  on  the 
entry  surface  of  the  target.  Material  strain-rate  depend¬ 
ence  is  considered  in  obtaining  the  associated  stresses  and 
the  work  rate  terms.  Certain  parameters  that  characterize  the 
field  of  stage  1  are  determined  by  utilizing  a  modification 
of  the  upper  bound  theorem  of  plasticity  to  include  dynamic 
effects. 

A  plastic  flow  velocity  field  is  also  assumed  for  each 
of  the  later  stages  with  parameters  introduced  to  describe 
the  geometrical  features  of  the  deformation  mechanism  under 
consideration.  In  these  cases,  the  parameters  are  obtained 
directly  from  the  mechanics  of  the  situation  without  the  need 
for  prior  test  information  or  minimization  procedure.  The 
resultant  flow  field,  in  turn,  determines  the  forces  that 
resist  the  perforation  process.  The  improved  engineering 
model  presented  in  this  paper  is  essentially  self-contained 
and  can  provide  adequate  analytical  treatment  of  the  various 
physical  phenomena  associated  with  perforation. 


Input  data  consists  of  the  projectile  diameter,  length, 
density,  and  impact  velocity,  and  the  mechanical  properties, 
density  and  thickness  of  the  target  material.  The  predic¬ 
tions  based  on  this  information  are  the  time  histories  of 
the  projectile  velocity  and  displacement  and  of  the  resisting 
force,  and  also  the  history  of  the  plastic  deformation  state 
of  the  target  plate  including  the  entry  "lip",  the  rear 
surface  bulge,  and  the  plug  shape.  In  addition,  the  analysis 
gives  the  exit  velocities  of  the  projectile  and  the  plug. 
Since  the  projectile  is  assumed  to  be  rigid,  flattening  is 
not  taken  into  account  although  the  general  method  can  be 
extended  to  include  this  effect.  A  flat  ended  projectile 
is  considered  for  simplicity,  while  other  nose  geometries 
would  require  modification  of  the  plastic  velocity  field 
but  would  not  cause  any  essential  difficulty  in  the  calcu¬ 
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General  Formulation 

The  general  procedure  for  analyzing  the  various  stages  of 
perforation  is  to  assume  a  plastic  flow  velocity  field  that 
describes  the  particular  mechanism.  This  field  would  be 
required  to  meet  the  kinematical  constraints  and  the  compati¬ 
bility  and  incompressibility  conditions  and,  in  some  cases, 
would  initially  contain  certain  undetermined  parameters. 

Those  parameter  values  would  be  obtained  by  minimizing  the 
expression  for  the  rate  of  plastic  work,  i.e.  the  energy 
dissipation  rate,  at  each  increment  of  penetration.  Convective 
inertia  terms  arising  from  the  motion  of  the  target  material, 
i.e.  terms  proportional  to  the  square  of  the  velocity,  are 
taken  to  be  part  of  the  dissipation  function  to  be  minimized. 

On  the  other  hand,  the  local  inertia  terms,  which  are 
proportional  to  the  acceleration,  are  treated  as  an  "effective 
mass"  of  material  acting  together  with  the  projectile  and 
therefore  not  subject  to  the  minimization  process. 

This  procedure  can  be  rationalized  on  the  basis  of 
Martin's  1972  theorem  on  acceleration  fields  [33  which 
states  that  the  rate  of  work  of  the  inertial  force  based  on 
an  assumed  acceleration  field  would  be  an  upper  bound  on 
the  actual  one.  Alternatively,  the  upper  bound  theorem  of 
plasticity,  in  the  absence  of  applied  forces,  can  be  inter¬ 
preted  as  applying  to  the  inertial  force  consisting  of  the 


projectile  and  "effective  mass"  of  target  material  multiplied 
by  their  respective  accelerations.  In  problems  involving 
prescribed  forces  acting  on  plastic  media,  e.g.  the  dynamic 
punch  problem,  the  material  acceleration  terms  would  be  part 
of  the  expression  to  be  minimized  based  on  direct  application 
of  the  upper  bound  theorem.  That  was  the  method  employed  by 
Tirosh  and  Kobayashi  in  [4]  . 

To  obtain  the  strain  rates  e ^ •  from  the  velocity  field, 
the  conventional  relation  is  used, 

Eij  -  (1/2)  <v.(j  ♦  vjf .),  where  lkk  =  0  (1) 

The  target  material  is  assumed  to  be  rigid-plastic  with  the 
flow  stress  a  function  of  strain  rate.  Work  hardening  can  be 
considered  indirectly  by  taking  the  flow  stress  to  be  the 
average  value,  with  respect  to  plastic  work,  over  the  full 
strain  range.  Strain  rate  dependence  of  the  flow  stress  is 
introduced  by  a  logarithmic  relation  which  has  been  used  by 
a  number  of  investigators,  e.g.  Lindholm  (51,  and  is 
convenient  for  the  perforation  problem. 


=  o  ( 1 +C  loa 
0  '10 


eeffl  = 


In  this  expression,  oq  is  the  reference  "static"  flow  stress 

•  Af  f 

(actually  corresponding  to  e  =  1) ,  the  material  constant 
C  is  taken  in  practice  to  be  zero  for  strain  rates  lower 


than  unity, 


and  a  is 

y 


l 
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the  flow  stress  corresponding  to  the  effective  plastic  strain 
rate  (calculated  in  the  usual  manner). 

Total  elastic  work  rate  in  a  volume  zone  V  can  be 

n 

expressed  as 


wf  =  /  a.  . dV  =  (a  )  £effV 
V  v  13  i]  y  n  n  n 


where  (a  )n  is  the  flow  stress  corresponding  to  the  average 

•ef  f 

effective  plastic  strain  rate  over  the  volume  V  ,  e  ,  defined 
r  n  n 


-  u/573>/v„l£  £~i 


.  .  £  .  .  dV 

lO  13 


For  zone  boundaries  with  a  tangential  velocity 
discontinuity  Av^  across  a  surface  Am  ,  i.e.  a  shear  velocity, 
the  plastic  work  rate  is 


When  the  surfaces  are  the  interfaces  between  flow  zones,  the 

stress  (o  )  is  taken  to  be  the  minimum  value  for  the  zones, 
y  m 

Friction  losses  due  to  relative  motion  between  the 
projectile  and  the  adjacent  target  material  can  also  be 
represented  by  (5)  with  a  friction  coefficient  u  that  could 
vary  between  0  and  1.  For  greater  generality,  these 
coefficients  are  taken  to  be  different  on  the  frontal  and 
lateral  faces  of  the  projectile  and  designated  and 
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respectively.  In  practice,  the  overall  influence  of  the 
frictional  forces  is  small  during  stages  1  to  4 .  They  are 
the  dominant  factor  for  stage  5  and  the  coefficients  are 
obtained  from  the  available  information,  e.g.  dynamic 
punching  tests. 

The  dissipative  work  rate  due  to  the  internal  stresses 

and  frictional  forces  would  then  be  quantities  W?l,  W™  and  W™ 

V'  s  f 

summed  over  the  volumes  and  areas  over  which  they  are  acting. 
These  quantities  can  be  non-dimensionalized  by  dividing  by 
irRsvoo  where  R  is  the  radius  of  the  projectile  (Fig. 2),  and 
v  is  its  current  velocity,  e.g. 

w‘  =  Ew"/(itR JV0  )  ;  w’  =  Ew“/(irR*vo  1  ;  etc.  (6) 

vnv  °  s  m  s  o 


The  work  rate  of  the  inertial  forces  in  the  zone  V 

n 

due  to  the  motion  of  target  material  is  given  by 


/  v.v.dV 

vn  1  1 


(7) 


where  w£  is  the  kinetic  work  rate  and  w”  is  the  dynamic 
work  rate  and  p  is  the  mass  density  of  the  target.  The 
total  acceleration  consists  of  the  local  and  convective 
components, 


°i  =  £(vi>  =  A(vi)  +  (vi'V)vi  (8) 

Upon  determining  v^  from  (8)  and  expressing  the  material 
velocities  in  the  plastic  zones,  v^,  as  functions  of  the 
projectile  velocity  v  and  the  geometry,  the  inertial  work  rate 
can  be  decomposed  in  the  form  of  (7) .  Non-dimensionalizing 
the  inertial  work  rate  terms  by  dividing  by  7rR*va0,  leads  to 


s 


Wk  =  (°tv2/oo)  (3] 


(9) 


w'd  =  (0tvR/oQ)  g 


(10) 


It  Is  noted  that  the  non-dimensional  kinetic  work  rate  is 

proportional  to  vJ  and  the  dynamic  work  rate  to  v.  The 

functions  g  and  g  in  (9,10)  depend  upon  non-dimensional 
12 

geometrical  parameters  of  the  flow  field. 

For  a  flat  nosed,  cylindrical  projectile,  the  work  rate 
of  the  total  force  acting  on  it  in  the  negative  velocity 
direction  is  given  by 


W  =  -  p  tt  R  L  w 
P  P 


(11) 


or,  in  non-dimensional  form. 


Wp  =  -  Pp  lv/po 


(12) 


where  pp  is  its  mass  density.  The  energy  rate  balance 
equation  for  the  complete  system  is  therefore 


•  I  »  I  i  I  •  I  *1  •« 

W  =  W,T  +  W  +  w.  +  w.  +  w. 

p  V  s  f  k  d 


or,  using  (10)  and  (12), 

-(ppL  +  PtR  g2Hv/oQ)  = 


where 


WT  =  WV  +  Ws  +  Wf  +  Wk 


(13) 

(14) 

(14a) 
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Equation  (14)  can  be  interpreted  as  a  simple  dynamical 
equation  of  motion  in  which  the  total  mass  is  that  of  the 
projectile  plus  an  effective  mass  of  target  material. 

Another,  and  different,  effective  mass  term  appears  in  the 
expression  for  .  In  this  respect,  the  analysis  is 
similar  to  that  of  (1J  but  the  various  terms  are  obtained 
from  more  fundamental  considerations. 

As  discussed  previously,  it  may  be  necessary  in  certain 
cases,  e.g.  stage  1,  to  obtain  certain  geometrical  parameters 
by  minimizing  the  energy  dissipation  terms  on  the  right 
hand  side  of  eq.  (14)  ,  W^.,  at  each  displacement  increment. 

This  will  determine  v,  and  a  simple  iterative  procedure  is 
used  to  obtain  the  projectile  velocity  v  and  time  for  the 
next  increment  in  displacement  Ax.  For  the  other  stages,  all 
the  geometrical  parameters  of  the  assumed  flow  field  are 
obtained  from  the  kinematic,  continuity  and  incompressibility 
conditions.  A  typical  complete  calculation  of  the  perforation 
process  requires  about  10  minutes  on  a  minicomputer, 
e.g.  Digital  PDP  11/44.  Detailed  descriptions  of  the  velocity 
fields  and  analytical  procedures  for  each  of  the  stages  are 
given  in  the  following  sections. 


Stage  1  -  Dynamic  Plastic  Penetration 


This  stage  would  develop  very  shortly  after  initial 
impact  and  is  characterized  by  continuously  growing 
plastic  zones  surrounding  the  projectile.  An  approximate, 
axisymmetric  velocity  field  that  satisfies  the  geometrical 
compatibility  conditions  was  suggested  by  Iddan  and  Tirosh 
[6]  for  the  punch  problem  and  is  used  in  this  analysis, 

Fig. 2.  The  velocity  components  in  cylindrical  coordinates 
for  the  various  flow  zones  and  the  corresponding  strain  rates 
from  eq.(1)  are  given  in  Table  1.  Velocities  normal  to  a 
zone  boundary  are  required  to  be  continuous  across  the  inter¬ 
face  while  discontinuities  in  the  tangential  direction  lead  to 
shear  energy  dissipation,  eq. (5) .  To  maintain  zero  volume 
change,  material  in  zone  III  must  move  opposite  to  the 
direction  of  penetration  leading  to  "frontal  petalling"  of  the 
entry  surface,  H^. 

The  longitudinal  extent  of  the  plastic  zone  ahead  of  the 
projectile  at  a  penetration  distance  x  is  given  by  aR  while 
the  plastic  zone  radius  is  r\R»  Fig. 2,  where  a  and  r\  are 
obtained  by  minimizing  W^,,  eq.  (14a),  at  each  increment  of 
penetration  Ax.  This  stage  will  end  when  the  plastic  zone 
depth  reaches  the  rear  surface  of  the  target  plate,  i.e. 

x  +  aR  =  H,  for  which  x«x  .  At  this  condition,  n  =  n  where 

o  o 

nQR  is  the  maximum  radius  of  the  plastic  zone  and  also  the  maximum 
extent  of  the  shear  zone  in  subsequent  stages,  so  that  the  shear 
zone  width  e  defined  in  (1,7]  would  correspond  to  R ( n  —  1 >  in  this 


treatment. 
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In  order  to  use  the  strain  rate  dependence  of  the  yield 
stress,  eq.  (2) ,  average  effective  strain  rates  are  calculated 
from  eq. (4)  for  zones  I  and  II  giving,  respectively. 


£eff  =  (v/aR)  (15) 

1 

f  1  ( n ) 

•  pf  f 

e  1  =  (v/aR)  -  (16) 

2  / 3(n2-D 1 


where 

f  (n)  =  2n2  -  /  n 4  +  3  ♦  n3  *n[  (1/3)  (/n^o  +  n3)]  (17) 

i 

For  zones  III  and  IV,  the  average  effective  strain  rate  is 
zero. 


Calculation  of  the  non-dimensional  energy  dissipation 
rate  terms  of  (14a)  for  the  velocity  field  described  in 
Table  1  leads  to  the  following  expressions: 

W  =  S  *  S  -  (18) 

v  1  3  /3(n’-i) 


where  f  (h)  is  given  by  (17),  and  S  and  S  are  the  strain 
1  12 

rate  factors  defined  by  (2)  for  zones  I  and  II  having  the 
strain  rates  (15)  and  (16)  respectively; 


•  I 

w_  = 


s  /3 (n2 -1 ) 


r  (4n*-5n2+l)  2x 

a[1+n  min(S  ,S  ))  +  -  +  —  V  (19) 

1  2  3an  R 
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and 


•  I 

Wf 


s  1  Uf 

3  /3  a 


2^s  n* 

/3 (n2-l ) 


(x+HL> 

R 


L 

6  +  -  (1-6) 
R 


(20) 


where  L  is  the  projectile  length  and  6=o  for  x£L-HL  while  6=1 
for  x<L-Hr .  The  frontal  lip  height  HT  is  the  accumulation  of 
the  incremental  lifts  of  material  in  zone  III  at  the  forward 
face  of  the  target,  where  AHL  =  Ax/(n2-1),  (see  Fig. 2). 


Expressions  for  Wk  and  Wd  are  given  by  (9)  and  (10),  and 

the  functions  g  and  g  for  the  assumed  velocity  field  are 

1  2 


g  (a, n)  =  - 

1  (n2 


Hn(n)  +  3-5 


g  (a, n)  =  - 

2  n2-i 


n2  f  1  f"8n2  (n2- 

-1 )  2  {  16a2  L  n2-1 

x  if  2n2  Jin  (n)  "1 

-  +  - - 1 

n2-i )  4a  I  nz-i 


„.l  .  i  ] 

3 


a 

_  + 
3 


(21) 


(22) 


Minimization  of  the  sum  of  the  dissipation  rate  terms 
W^, ,  (14a),  is  perurmed  at  each  increment  of  penetration  Ax 

with  respect  to  the  parameters  that  describe  the  plastic  zone, 

•  * 

a  and  n.  A  workable  procedure  is  to  formally  minimize  WT  with 
respect  to  a  which  leads  to  a  cubic  in  a  with  coefficients  that 
are  functions  of  n  and  S n  ( e ) .  A  numerical  search  technique  is 
then  employed  to  find  values  of  a  and  n  that  would  satisfy  the 
minimization  conditions. 
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Stage  2  -  Bulge  Formation 

Once  the  front  of  the  plastic  zone  developed  in  stage  1 
reaches  the  rear  free  surface  of  the  target  plate,  x+aR=H 
for  x=xq,  the  mode  of  plastic  deformation  transforms  into 
one  characterized  by  the  formation  of  a  spherical  bulge  on 
that  surface.  Uplifting  of  material  at  the  impact  surface 
then  ceases  and  the  material  in  a  truncated  conical  sector 
forward  of  the  projectile  is  pressed  as  a  spherical 
expansion  to  form  the  bulge.  This  effect  is  clearly  observed 
in  ballistic  tests,  e.g.[7).  The  essential  physics  of  this 
mode  is  that  of  target  mass  movement  compatible  with 
incompressibility.  The  associated  plastic  flow  velocity 
field  determines  the  resisting  forces  acting  on  the 
projectile.  No  minimization  procedures  are  therefore 
required  for  this  stage. 

In  order  to  analyze  the  velocities  and  forces,  it  is 
convenient  to  introduce  spherical  coordinates,  r,  0,  <)>,  with 
origin  0  determined  by  the  semi-cone  angle  8  for  the  region 
forward  of  the  projectile,  i.e.  zone  II  of  the  proposed  flow 
field  shown  in  Fig. 3.  A  similar  radial  expansion  flow  field 
has  been  used  by  Avitzur  [8]  and  others  in  analyses  of  metal 
processing.  The  relation  between  the  spherical  bulge  geometry, 
and  the  movement  of  the  projectile  is  obtained  from  the 
incompressibility  requirement,  narn*ly. 


i 


»n3R3  llzs.ose)  7  ,(2tcosg)  =  ^ 

b  3sin 3  8  0 


where  nbR  is  the  radial  extent  of  the  bulge  in  the  plane  of 

the  rear  surface  of  the  plate,  and  x-x^  is  the  displacement 

of  the  projectile  from  the  end  of  stage  1  (Fig. 3).  A  further 

geometrical  relationship  can  be  obtained  between  n.  and  8  from 

b 

Fig. 3, 


nb  =  t (H~x) / R]  (tanB) +  1 


which,  combined  with  (23),  leads  to  an  equation  relating  8 
to  (x-xQ), 


(H~x)  (tanB)  +  1 


3 (sin’ 8) (x-*o) 
R(1-cos8)?  (2+cosB) 


This  equation  can  be  solved  numerically  for  8  for  a  given 
projectile  displacement  x .  The  bulge  height  in  the  axial 
direction,  H,  ,  can  also  be  obtained  from  the  geometry. 


H,  =  (n.R/sinB)  (1-cosB) 
b  b 


In  zone  II,  the  only  non-zero  velocity  is  vf  in  the 
radial  direction  r. 


vr  =  v ( r  /r ) 2  cost 


(27) 


where  v  is  the  projectile  velocity  and  r  and  r ^  are  indicated 
in  Fig. 3.  For  this  velocity  field,  the  strain  rate  components 
from  ( 1 ) ,  are 


=  -2v (r2/r 3 )  cosP 

(28) 

1 

=  c,,  =  v(r2/r3)  cos0 

<p<p  i 

(29) 

=  -(1/2)v(r2/r3)  sinB 

(30) 

Based  on  the  velocities  and  strain  rates  in  zone  II,  the  work 
rates  Wy,  Wd  and  could  be  calculated  for  that  zone  as  they 
were  obtained  in  stage  1  using  eqs  .( 3 , 4 , 7 , 8 )  . 

Zone  I  is  considered  to  be  a  "dead"  zone  with  the 
material  moving  uniformly  with  the  projectile  velocity  v  so 
that  a  tangential  velocity  discontinuity  exists  on  the 
interface  between  zones  I  and  II.  The  shear  work 
rate  W^12^  due  to  the  velocity  jump  across  the  boundary 
between  zones  I  and  II  is  calculated  from  eq. (5) .  Zone  I 
also  contributes  a  work  rate  but  not  and  Wv  since  the 
velocity  is  constant  throughout  that  volume. 

Zone  VI  is  unaffected  by  the  flow  process.  The  boundary 
of  III,  IV  and  V,  nQR,  is  the  value  for  nR  obtained  at 
the  end  of  stage  1,  when  x=xQ.  The  velocity  field  for  zone 
V  is  not  specified  since  its  influence  on  the  stage  2 
process  would  be  negligible  except,  possibly,  for  frictional 
forces  acting  along  the  projectile  surface.  The  work  rate 
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of  the  frictional  forces  Wf  is  again  calculated  by  the  methods 
described  for  stage  1 .  Zone  IV  is  a  narrow 
region  at  the  rear  surface  and  its  contributions  are 
negligible. 

A  proper  velocity  field  for  zone  III  of  Fig. 3  would  be 
difficult  to  construct.  Instead,  the  overall  influence  of 
this  region  is  considered  in  an  approximate  manner.  From 
observations  of  cavity  cross-sections,  e.g.  [7,  9],  it  is  seen 
that  this  region  acts  primarily  as  a  transition  zone  which 
undergoes  shearing  deformation  due  to  the  material  velocity 
along  the  interface  with  zone  II  (assuming  no  discontinuity 
in  that  velocity)  and  the  absence  of  material  movement  in 
zone  VI . 

The  shear  strain  rate  in  zone  III,  y  ,  is  assumed  to  be 

3 

constant  along  a  circular  arc  d(r').  Fig. 3,  and  is  approximated 
as  the  velocity  at  the  interface  with  zone  II,  v  ,  divided  by 

2  3 

the  radial  width  of  the  zone  described  by  the  arc  length  d  so 
that 

=  (Rzv  cot6)/[(r,)zB(n0R-r*sin6))  (31) 

where  r'  is  measured  from  the  origin  0  along  the  interface 
surface  A  .  On  the  basis  that  the  plastic  work  rate  in  the 

2  3 

«3ne  III  volume  is  due  entirely  to  shearing  deformation,  the 

term  can  be  expressed  by  eq.  (3)  where  =  (y  )_.//!. 

V  33  avg 

The  average  shear  strain  rate  in  zone  III  can  be  calculated 
by  integrating  (31)  over  the  volume  of  zone  III  and  dividing 
by  the  volume  V  . 
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Since  the  detailed  velocity  field  of  zone  III  is  not 
specified,  the  work  rates  due  to  the  dynamic  effects,  and 
W^,  are  approximated  by  the  work  rate  of  equivalent  body 
forces  (per  unit  volume)  acting  on  the  target  material  in 
zone  III.  Those  body  forces  are  assumed  to  be  proportional 
to  y  as  an  approximation  to  both  the  local  and  convective 
acceleration  effects.  The  coefficient  of  proportionality, 

K^,  is  obtained  by  enforcing  continuity  of  the  total  resisting 
force  acting  on  the  projectile  at  the  time  corresponding 
to  the  transition  from  stage  1  to  2.  Although  has  the 
dimension  of  dynamic  viscosity,  it  is  not  a  material  property 
but  a  parameter  that  represents  overall  inertial  resistance 
to  flow  in  an  approximate  manner.  In  the  absence  of  other 
information,  is  taken  to  be  constant  throughout  stage  2 
and  the  subsequent  stage  3.  The  inertial  work  rate  term  for 
zone  III  is  therefore  approximated  as 

W(J)  *  /( K.y  )  y  dV  (32) 

I  V  b  3  3 

In  actual  computations  for  relatively  thick  target  plates, 
the  contribution  of  (32)  is  small  compared  to  that  for  Wy3^ 
so  the  approximate  nature  of  (32)  is  not  significant. 

The  total  work  rates  for  stage  2  would  therefore  consist 
of  the  dynamic  term  from  zone  I,  the  shear  term  from  the 
interface  of  zones  I  and  II,  the  various  terms  from  zone  II, 
the  shearing  and  inertial  terms  for  zone  III,  and  the 
frictional  term  from  zone  V.  These  would  then  be  used  in 
the  energy  balance  equation  (13)  to  obtain  the  motion  and 
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forces  at  each  increment  of  projectile  displacement.  For  this 
stage,  the  velocity  field  is  completely  specified  by  the 
requirements  of  incompressibility  and  geometry  so  that  the 
determination  of  parameters  by  a  minimization  procedure  is  not 
required.  Stage  2  continues  until  x=xb  when  the  bulge  width 
n,  R  reaches  a  maximum  value,  nPR.  Further  continuation  of  the 

□  D 

deformation  mode  would  lead  to  physically  unrealistic  situations 
which  indicates  that  another  deformation  mechanism  initiates  at 


that  condition. 
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Stage  3  -  Bulge  Advancement 

Guided  by  observations  of  bulge  formation  and  growth, 
e.g.  [7,9]  the  alteration  to  the  deformation  pattern  is 
taken  to  be  an  axial  displacement  of  the  bulge  region  with 
the  new  bulge  radius  n^R  ending  on  the  former  bu’ge  surface. 
Fig. 4.  This  bulge  advancement  is  due  to  target  material 
displaced  by  the  projectile,  and  the  resulting  geometry  is 
determined  by  the  incompressibility  condition. 

From  Fig. 4,  an  increment  of  projectile  displacement 
Ax  leads  to  an  advance  in  the  bulge  where  the  new  spherical 
increment  has  a  sector  angle  3  and  the  coordinates  of  its 
intersection  with  the  former  bulge  surface  are  n^R  and  £. 

The  relation  between  the  geometrical  parameters  is 

(1-cosB)2  (2+cosB)  (1-cosa)2  (2+cosa)  3Ax 

- - - - - —  =  o  (34 

sin  8  sin  a  n£R 

where  a  is  the  radial  angle  from  the  former  sector  origin 
0“  to  the  new  boundary  limit  n^R,  and  is  given  by 


a 


sin-1  I(rib/n^)  sinB” ] 


(35) 
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and  nb  and  6  are  geometrical  parameters  from  the  previous 
step.  The  other  necessary  geometrical  relations  are 

B  =  tan-1 [R(nb-1)/ (H+C-x)]  (36) 

6  =  C  +  (n^R/sinB  ) (cosa-cosB  )  (37) 

and  the  current  bulge  height  is  given  by 

Hb  =  £  *  (nbR/sinB) (1-cosB)  (38) 

These  equations  are  sufficient  for  obtaining  the  bulge 
geometry  at  each  increment  of  projectile  displacement.  The 
work  rates  of  the  various  forces  associated  with  the  plastic 
flow  field  could  be  calculated  by  the  procedures  described 
for  stage  2  and  these  stages  will  be  fully  continuous. 

Stage  3  of  bulge  advancement  could  continue  even  for  x>H. 


Failure  Mechanisms 


The  penetration  processes  described  by  stages  1  to  3  are 
common  to  a  large  range  of  ballistic  impact  conditions,  but  the 
terminal  stage  leading  to  failure  of  the  target  and  exit  of  the 
projectile  is  characterized  by  a  number  of  possible  mechanisms. 
Some  of  these  are  described  and  illustrated  in  [9].  It  is 
possible  for  stages  2  or  3  to  lead  directly  to  a  terminal  stage. 
This  could  occur  by  the  development  of  adiabatic  shear  bands  on 
the  interface  A  of  zones  II  and  III,  (Figs. 3,4)  or  to  the 

2  3 

target  material  being  relatively  brittle  and  failing  in  shear 
on  that  surface.  Such  early  failures  in  the  penetration  process 
lead  to  separation  of  a  truncated  conical  sector  (a  conical  plug) 
consisting  of  zones  I  and  II,  (Figs. 3, 4)  from  the  target  plate 
and  are  "low  energy"  modes  of  plugging  failure  as  discussed, 
e.g.  in  [10] . 

Alternatively,  and  more  commonly  for  ductile  targets  that 
do  not  fail  by  adiabatic  shearing,  a  conventional  plugging 
mechanism  of  "intermediate"  energy  (with  a  cylindrical  plug) 
could  develop  from  stages  2  or  3  by  the  formation  of  a  flow 
field  with  a  shear  zone  within  the  maximum  bulge  radius.  For 
this  mode,  the  radius  of  the  ejected  plug  is  essentially  equal 
to  or  less  than  that  of  the  projectile. 

Very  ductile  target  plates,  which  exhibit  "high  energy" 
ballistic  resistance  could  experience  failure  by  mechanical 
instability  of  the  system  when  stage  3  is  carried  out  to  its 
geometrical  limitation.  The  resulting  plug  is  relatively 
short  and  has  a  spherical  cap  which  is  the  frontal  part  of 
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the  bulge.  This  and  the  preceding  modes  of  failure  are 
described  in  more  detail  in  the  following  sections.  They  are 
the  ones  most  commonly  observed  in  ballistic  tests  on  metallic 
plates,  but  other  failure  mechanisms  and  combinations  of  the 
above  modes  are  also  possible. 

a.  Adiabatic  Shearing 

The  appearance  of  adiabatic  shear  bands  on  surface  A 
during  stages  2  or  3  would  lead  to  a  velocity  discontinuity 
(slip)  along  A  (Figs. 3,4)  thereby  separating  the  interior  of 

2  3 

zone  III  from  the  resisting  force  field.  A  physical  descrip¬ 
tion  of  the  process  is  given  in  [10],  and  an  analytical 
treatment  of  the  formation  of  adiabatic  shear  bands  in  metals 
subjected  to  rapid  shear  deformations  is  presented  in  [11]. 

It  is  possible  to  modify  the  calculations  for  stages  2  and  3 
for  this  condition  so  that  a  flow  field  analysis  leading  to 
failure  can  be  obtained.  In  this  analysis,  the  plastic  and 
inertial  work  rates  of  zone  III  are  neglected  and  a  shear  work 
rate  on  A  is  considered.  From  (5)  and  (6) ,  this  term  is 

2  3 

W'(?3)  =  (2£nn. )/ (/3tan6)  (39) 

s  b 

Shear  strains  would  accumulate  in  the  band  region  during 

stage  3  and  an  approximate  measure  of  the  shear  strain  at  the 

rear  surface  is  y  :  5/ [R (n^-n^) ] .  Failure  will  occur  at 

x  =  x  when  y  =  y  where  y  ~  /5r  .  The  mass  of  the 

s  cr  cr  cr 

ejected  target  material  (plug)  consists  of  zones  I  and  II, 

Fig. 4,  at  the  failure  condition  and  its  average  velocity 
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could  be  obtained  by  momentum  considerations  from  the  velocity 
field.  Initially,  the  projectile  and  plug  will  move  together 
but  they  will  separate  as  the  projectile  is  slowed  due  to 
frictional  forces. 

b.  Brittle  Failure  in  Shear 

Shear  failure  along  the  interface  could  take  place 

when  the  material  rupture  condition  is  reached  along  that 
surface  during  stages  2  or  3.  There  is  no  velocity  disconti¬ 
nuity  along  A ^  prior  to  failure  by  this  mode  so  the 
standard  methods  for  stages  2  and  3  are  used  and  an  appropriate 
rate  dependent  failure  criterion  is  introduced  which  involves 

integrating  (>■)_.„  over  the  time.  Upon  the  critical  condition 
3  avg 

being  realized,  a  truncated  conical  plug  will  eject  as  in  the 
previous  case. 

c.  Plug  Formation  and  Exit  (Stages  4  and  5) 

Conventional  plug  formation  could  commence  at  any  time 
during  stages  2  and  3  and  has  to  be  examined  in  parallel  with 
the  bulging  processes.  The  mode  having  the  lower  resistance 
to  continued  penetration  will  be  the  active  one.  Plugging 
could  occur  in  a  variety  of  geometries  and  could  also  combine 
with  scabbing  to  form  complex  failure  modes  as  indicated  by 
Fig.3e  of  [91.  A  complete  analysis  would  require  consideration 
of  each  geometry  but  only  the  simple  cylindrical  plug  is 
described  in  this  section. 

The  flow  field  assumed  for  plug  formation  (stage  4)  is 
shown  in  Fig. 5  alongside  the  field  for  stage  3.  For  this  mode, 
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all  the  target  material  in  zone  I  forward  of  the  projectile 

is  assumed  to  move  at  the  projectile  velocity  v  and  is 

restrained  by  shear  stresses  in  zone  II.  Zone  III  is 

considered  to  contain  fractured  material  and  to  be  ineffective 

in  shear.  The  coordinate  C  defining  the  boundary  between 

zones  II  and  III  is  obtained  by  setting  the  shear  strain  at 

that  point  due  to  the  prior  shear  deformation  in  stage  3, 

approximated  by  c/[R(no-1)l,  equal  to  the  ultimate  strain 

Y  .  Those  zones  are  assumed  to  be  limited  in  width  to  the 
'cr 

maximum  bulge  diameter  n^R  on  the  basis  of  physical 
observations  so  that  the  shear  strain  rate  along  the  inter¬ 
face  of  tb**  plug  and  the  restraining  zone  is 

y  z  v/[R(n£-l>]  (41) 

If  transition  to  plugging  occurs  during  stage  2,  then  the 
current  value  of  the  bulge  width,  nfaR  is  used  instead  of 
n°R  in  (41) . 

Once  the  plugging  mode  has  formed,  further  projectile 
advancement  is  resisted  by  the  work  rates  of  the  inertial 
forces  of  zones  I,  II  and  III,  the  shear  work  rate  in  zone  II, 
and  by  friction  along  the  lateral  surface  of  the  projectile. 
These  quantities  are  determined  in  a  manner  similar  to  those 
of  stages  2  and  3.  In  calculating  the  work  rates  of  the 
inertial  forces  in  zones  II  and  III  (Fig. 5)  the  parameter 
that  represents  dynamic  effects  in  stages  2  and  3,  [eq. (32) ) , 
must  be  modified  to  account  for  the  difference  in  the  velocity 
field.  A  similarity  argument  is  used  to  obtain  a  new  value  of 


I 


i 
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K^for  this  stage. 


The  plug  length  H 


P 


is  obtained  from  the  geometry. 


H  =  R(cot(j>  -  cote  )  (42) 

P  P 


where 


=  sin-1(sin8  /n.  ) 
P  b 


(43) 


and  6^  is  the  sector  angle  for  the  flow  field  of  stage  2  or  3 
at  the  onset  of  plug  formation,  x  =  x  . 


The  bulge  will  continue  to  advance  until  the  ultimate 
shear  strain  is  realized  in  20ne  II.  At  this  condition, 
x  =  x  where  x  =  x  +  H  and  H  =  Y_  RtnP-l).  The  plug  will 

S  S  p  6  0  Ci  0 

then  detach  and  exit  with  the  projectile.  A  distinct  exit 
stage  (stage  5)  can  be  considered  consisting  of  the  projectile 
subjected  to  frictional  forces,  and  a  free  cylindrical  plug 
with  a  cap  corresponding  to  the  bulge  shape. 


d.  Ductile  Failure 

If  failure  by  any  of  the  previously  described  or  alternative 
modes  does  not  occur,  then  stage  3  would  continue  until  a 
mechanical  instability  condition  is  reached.  This  arises  when 
the  projectile  displacement  beyond  the  rear  surface  of  the 
target  exceeds  the  bulge  advancement  coordinate  £  leading  to 
high  shear  stresses  on  the  surface  A^,  Fig. 4.  It  is 
reasonable  to  expect  that  failure  would  then  precipitate  on 
and  that  a  plug  will  form  consisting  of  zones  I  and  II  at  that 
condition  (x>H)  which  will  have  a  spherical  cap  slightly 


wider  than  the  projectile.  This  capped  plug  would  stay 
attached  to  the  projectile  as  it  exits  from  the  plate  and 
will  subsequently  separate. 
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Numerical  Examples 

Numerical  exercises  based  on  variations  of  the  ductile 
plugging  exit  mode  were  carried  out  for  the  cases  of  A.P. 
projectiles  ot  7.62  mm  diameter  reported  in  [7]  and  [12]. 
Ductile  plugging  was  observed  in  those  tests  to  be  the 
dominant  exit  mechanism  and  was  accompanied  in  a  few  instances 
by  scabbing  and  fragmentation.  Only  the  hard  core  of  the  A.P. 
projectile  was  considered  for  the  calculations.  The  core  was 
idealized  to  be  a  rigid  cylindrical  rod  having  the 
maximum  core  diameter  of  6.1  mm  and  the  actual  core  length  of 
23.8  mm.  The  strain  rate  parameter  in  (2)  was  taken  to  be 
C  =  0.025  based  on  tensile  test  results,  and  the  friction 
coefficients  were  set  to  be  =  0.10  and  =  0.05.  A  value 
of  p  =  0.10  is  given  in  the  literature  for  dynamic  metal 
working  conditions  and  a  lower  value  should  be  used  here  for 
the  lateral  surfaces  because  of  the  higher  velocities  and 
temperatures. 

The  computed  and  test  results  are  listed  in  Table  2, 
from  which  it  appears  that  the  predictive  capability  of  the 
present  analysis  is  fairly  good.  Some  further  explanations 
are  required  for  reading  Table  2.  All  the  test  results 
listed  are  given  in  [12]  and  those  of  tests  1,  2,  3  are  also 
presented  in  [7] .  For  some  reason,  the  value  of  b  for  test 
3  in  [7]  is  not  consistent  with  the  correct  value  given  in 
[12]  and  in  Table  2.  The  width  of  the  shear  zone  in  the 
tests,  e,  was  measured  to  be  the  maximum  extent  of  the  visible 
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shear  strains,  e.g.  Figs.  5,  7  of  [7],  and  the  corresponding 
quantity  in  the  present  analysis  is  R(nQ-1)  listed  in  the 
last  column  of  Table  2.  Examination  of  cavity  cross-sections, 
e.g.  Figs.  5,  7  of  [7],  were  used  to  obtain  b,  the  length  of 
the  shear  zone  within  the  original  target  thickness.  The 
corresponding  parameter  in  the  analysis  is  H-xq  where  xq  is 
the  extent  of  the  penetration  at  the  end  of  stage  1.  Predicted 
plug  lengths,  H^,  obtained  from  the  computations  for  stage  4 
are  listed,  but  test  data  for  the  actual  plug  lengths  is  not 

available.  As  expected,  <  (H-xq)  in  this  analysis  whereas 

b  is  assumed  in  (i]  to  correspond  to  the  plug  length. 

Test  results  for  the  final  velocity  of  the  projectile 
give  a  range  of  values  which  are  a  consequence  of  the  ejected 
plug  and  fragments  activating  the  measurement  system  prior  to 
the  projectile  in  some  of  the  tests.  Calculated  values  for 
the  final  projectile  velocity  based  on  the  present  analysis 
are  generally  within  the  experimental  range  and  appear  to  be 
more  consistent  with  the  test  results  than  those  of  [1,7]. 
Predicted  plug  velocities  are  slightly  higher  than  those  of 

the  projectile  because  of  the  friction  on  its  lateral  surface 

during  stage  5. 

In  addition  to  the  primary  exit  modes  discussed  previously, 
there  are  secondary  differences  in  the  exit  mechanisms.  A 
number  of  possible  variations  of  the  basic  ductile  plugging 
mode  have  been  incorporated  in  the  calculation  procedure  and 
the  one  appropriate  for  the  test  was  used.  It  is  recognized 
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that  this  procedure  is  an  empirical  input  and  corresponds,  in 
effect,  to  distinguishing  fine  differences  in  the  exit 
behavior.  These  differences  appear  to  be  related  to  the 
observed  scatter  of  terminal  velocity  values  for  nominally 
identical  test  conditions  which  is  especially  pronounced  near 
the  ballistic  limit.  The  variations  of  the  ductile  plugging 
mode  that  were  considered  are  as  follows: 

PQ  -  a  fully  ogival  shaped  plug  within  the  cylindrical 
boundary. 

P  -  a  cylindrical  plug  with  an  ogival  nose,  e.g.  Fig.  3b 
of  [91. 

P  -  a  fully  cylindrical  plug  as  described  in  stage  4  of 

2 

the  analysis,  e.g.  Figs.  4,  8  of  [7],  Fig.  3a  of  [9], 

P  -  a  cylindrical  plug  with  a  wider  head  (combined 

plugging  and  scabbing)  due  to  the  combination  of 
zones  I  and  III  in  stage  4  (Fig.  5),  e.g.  Fig  3e 
of  [9]. 

Two  plugging  conditions  were  calculated  for  test  7  since  the 

actual  shape  seemed  to  have  features  of  both  P  and  P  .  The 

12 

results  are  an  example  of  the  range  in  final  velocity  values 
obtainable  from  variations  in  the  basic  plugging  mode. 

Calculated  histories  of  the  projectile  displacement, 
velocity,  and  resisting  force  are  shown  in  Figs.  6  and  7  for 
tests  2  and  3.  For  the  thinner  target,  test  2  -  Fig.  6,  the 
duration  of  the  high  resisting  force  during  stage  1  is  very 
short,  1.6  us,  and  stages  2  and  3  play  an  important  role.  A 
rapid  drop  in  the  resisting  force  occurs  during  stage  4  (at 
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t  =  10.3  us,  x  =  7.4  mm)  due  to  separation  of  the  assumed 
cylindrical  plug  ( P ^ )  from  the  target  plate  at  that  instant. 
Stage  4  continues  until  the  plug  is  fully  ejected  from  the 
target  plate.  The  maximum  value  of  x,  33.4  mm,  is  the 
distance  for  complete  exit  of  the  projectile  from  the  deformed 
target  at  the  end  of  stage  5.  Calculated  items  not  listed  in 
Table  2  for  test  2  are:  plug  mass  0.6  g,  maximum  bulge  height 
3.7  mm,  maximum  bulge  diameter  10.1  mm.  Effective  strain 
rates  in  the  bulge  region  and  along  the  surface  A^  ^  of  stages 
2  and  3  were  of  the  order  of  10s  sec-1.  Test  3,  Fig.  7,  was 
for  a  thicker  target  plate  and  the  duration  of  stage  1  is 
longer,  8.5  us.  Again,  stages  2  and  3  are  important  factors. 
Plug  separation  takes  place  at  t  =  22.1  ps  when  x  =  13.7  mm 
but,  because  of  the  assumed  ogival  plug  shape  (PQ) ,  there  is 
no  sharp  drop  in  the  resisting  force  during  stage  4.  Additional 
computed  quantities  for  this  case  are:  plug  mass  0.5  g,  maximum 
bulge  height  3.8  mm,  maximum  bulge  diameter  10.5  mm.  Effective 
strain  rates  were  about  1/3  lower  for  the  thicker  target  plate. 
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Particle  Velocities  and  Strain  Rates 
Flow  Field  Assumed  During  Stage  1  of 


for  Each  Zone  of  the  Plastic 
the  Penetration  Process 


Experimental  Results  and  Measurements  Analytical  Predictions 
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Pig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


Fig.  5 


Fig.  6 


Fig.  7 


Caption  for  Figures  . 

-  Schematic  of  five  stage  sequence  of  deformation 
mechanisms  for  perforation  process. 

-  Plastic  flow  field  for  stage  1  of  perforation 
process  -  dynamic  plastic  deformation. 

-  Plastic  flow  field  for  stage  2  of  perforation 
process  -  bulge  formation. 

-  Plastic  flow  field  for  stage  3  of  perforation 
process  -  bulge  advancement. 

-  Plastic  flow  field  for  stage  4  of  perforation 
process  -  plug  formation. 

-  Predicted  histories  of  projectile  displacement, 
velocity,  and  resisting  force  during  perforation 
process  (Test  No. 2,  Table  2). 

-  Predicted  histories  of  projectile  displacement, 
velocity,  and  resisting  force  during  perforation 
process  (Test  No. 3,  Table  2). 
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